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Abstract. We carried out a simulation of heavy ion collision using a time-dependent density 
functional theory. We call it the canonical-basis time-dependent Hartree-Fock-Bogoliubov 
theory (Cb-TDHFB) which can describe nuclear dynamics in three-dimensional coordinate 
space, treating nuclear pairing correlation. We simulate 20 O+ 20 O collision using the Cb-TDHFB 
• with a contact-type pairing functional, and show the behavior of gap energy which is decreasing 

' and vibrating while colliding. 

>. 

£SJ , 1- Introduction 

Atomic nucleus is a self-binding finite quantum many-body system whose size is several 
femtometers (= fm = 10 -15 m). It is composed of two kinds of fermions (protons and neutrons) 
with spin 1/2 that are called "nucleons". The number of nucleons in a nucleus is a few hundreds 
at most. There are about 3,000 nuclei observed in experiments, and yet more to be discovered. 
£Sj | In order to study the structure of nuclei, we investigate many excited modes of many nuclei. To 

study nuclear collision microscopically is also very important. Nucleus has "magic numbers", 
which are analogous to those in atoms. This suggests that the nucleon in the nucleus is freely 
moving in an attractive average potential. The nucleons create the mean field by themselves. 

In the time-dependent mean-field theory, the heavy-ion collision is expressed as colliding of 
two mean fields. We should investigate the behavior of single-particle states while colliding two 
fields, to expect the structure of nuclei after the heavy-ion collision. Nuclear pairing correlation 
has a significant role to define single-particle states. It allows the fractional occupation 
probability of single-particle states due to the pair scattering, which is important to reconstruct 
a new mean field. The time-dependent Hartree-Fock theory (TDHF) is well-known as a useful 
tool to study nuclear dynamics [lj, however it can not include effects of pairing correlation. 

There is a theory to treat pairing correlation in nuclear dynamics self-consistently, which is 
called the time-dependent Hartree-Fock-Bogoliubov theory (TDHFB). But, there is no study 
of heavy ion collision using TDHFB with a modern effective interaction in three-dimensional 
space, because the study is hard to prepare converged initial state, and it needs the huge 
computational resource [2j. We proposed a new TD mean- field theory to study nuclear dynamics 
treating pairing, which is named the canonical-basis TDHFB (Cb-TDHFB) [3]. The Cb-TDHFB 
is derived from full TDHFB equations represented in canonical basis which diagonalizes density 



matrix, with Bardeen-Cooper-Schrieffer (BCS)-like approximation for pairing functional (pp,hh- 
channel). We confirmed that linear response calculations of the Cb-TDHFB are a good 
approximation of the quasi-particle random phase approximation (QRPA) which is a small 
amplitude limit of TDHFB [3]. 

In the present work, we apply the Cb-TDHFB to heavy ion collision in the three-dimensional 
Cartesian coordinate space using a modern effective interaction (Skyrme: SkM*) with a contact 
pairing. In the followings, we introduce the Cb-TDHFB equations and the contact pairing 
functional. Then, we show results of 20 O+ 20 O which indicate an collision of the gap-energy. 



2. Formulation and Procedure 

2.1. Basic equations and contact pairing 

The Cb-TDHFB equations can be derived from the TDHFB equations represented in canonical 
basis with a simple approximation for pairing functional [5]. The TDHFB equations can be 
written in terms of the generalized density matrix TZ(t) and the generalized Hamiltonian 
1Z(t) composes the one-body density matrix p^t) = {^{t)\clc^\^{t)) and the pairing-tensor 
Kfj,v(t) = (^(t^Cj/Cftl^^)): fj,,v mean arbitrary basis. T-L(t) has h(t) and A(t) corresponding to 
single-particle Hamiltonian and pair potential, respectively. We can express TDHFB state at 
any time in the canonical (BCS) form as 

\m) = U(Mt)+vi(t)4(t)c\(t))\Q), (i) 

where ui(t),vi(t) are time-dependent BCS factors and cj is a creation operator of canonical basis 
which diagonalizes density matrix pit). Then, we choose the BCS form of pair potential as 

A,(t) = (2) 

fc>0 

where Kk(t) = Uk(t)vk(t) corresponds to the pair tensor n(t) in the canonical-basis and "P^ife is 
the anti-symmetric two-body matrix element. The subscripts I and k mean the pair of canonical 
basis I and k, respectively. We can obtain the Cb-TDHFB equations from the canonical form 
Eq-dU) of TDHFB, and with the pair potential Eq.@, as follows. 

(h(t) - Vl (t))Ut), ih?&& = (h(t) - m {t))Utl 

Ki(t)AUt) - Kt(t)Ai(t), 

( Vl (t) + m (t))Ki(t) + Mt)(2 Pl (t) - 1), (3) 

where rji(t) = (4>i(t) \ h(t) \ (j>i(t)) + ih(^^ \ 4>i(t)). The Cb-TDHFB equations compose the time- 
evolutions of the canonical basis (j)i(t), <j>i(t), the occupation probability pi(t) = \vi(t)\ 2 and the 
pair probability Ki(t). They conserve orthonormal property of the canonical basis and average 
particle number. When we choose a special gauge condition rji(t) = ej(t) = {<fii(t)\ h(t) \ 4>i)(t), 
they conserve average total energy. 

We introduce neutron-neutron and proton-proton BCS pairing. The BCS pairing matrix 
elements VT T , T is written as 

ll,kk 

V ll,kk = I dridr 2 ^ $ (r u ai)^f(r 2 , a 2 )V T (n, af, r 2 , a 2 ) 

0-1,0-2 

x [(j) k (r 1 ,a 1 )(f>) : {r 2 ,a 2 ) - <t>k(ri,a 1 )(/) k (r 2 ,a 2 )} . (4) 
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We introduce the spin-singlet contact interaction to Eq.(JU): 

V T (n,(Ti; r 2 , (7 2 ) = Vq~ — — S(rt - r 2 ), (5) 

where r indicates neutron or proton channel and Vq is a strength of pairing[5]. In this present 
work, we do not use the time-reversal relation between /- and /-canonical basis. 

2.2. Procedure and calculation space 

At first, we prepare the initial states (HF or HF+BCS ground states) with the wave functions 
for projectile and target at some impact parameter b at distance where they feel only Coulomb 
interaction. Then, we boost the wave functions and translational calculate the time-evolution of 
nuclear densities to obey Eq.([3]). In this present work, we use the three-dimensional Cartesian 
coordinate-space representation for the canonical states, cf>i(r, a; t) = (r, a\ <t>i{t)) with a = ±1/2. 
The coordinate space is a box whose x-, y-sides are 30 fm and z-side is 50 fm, discretized in the 
square mesh of Ax = Ay = Az = 1.0 fm. 

3. Results 

We simulate the 20 O+ 20 O collision with incident energy E\ n x = 40 MeV which is higher than 
the Coulomb barrier in this system. Initial distance between projectile and target is 12 fm 
with 6 = fm. 20 O (Z=8, N=12) has superfluidity for neutrons. The average gap-energy 
A = J2i>o ^i/J2i>o i s 1-901 MeV. The number of canonical basis is 24(1.2 times that of TDHF). 

Figure 1 shows the time evolution of neutron density distributions. After the touch 
(t = 98.66 fm/c: panel (b) in Fig.l), the density-distribution vibrates keeping the prolate shape. 
This shape of a new mean field corresponds to the shape of 40 S (Z=16, A^=24) ground state in 
intrinsic frame. Figure 2 shows the time-evolution of neutron gap-energy for each canonical basis 
defined in Eq.Q. From the touching configuration, the gap energies are decreasing significantly. 
The decrease is similar to the results of old work [6]. We can see the vibration of gap energies 
related with nuclear density behavior. 

We simulated the 20 O+ 20 O collision using time-dependent mean field theory including the 
effects of nuclear superfluidity, and showed the vibration of gap energies in the collision. 
Currently, we are investigating the detail of gap energy vibration and relation with the excited 
states of 40 S. To study the pairing effects for collision phenomena, we should improve our 
program code to parallel one. We parallelize the code using MPI for the canonical basis, because 
the canonical basis conserve orthogonality in the Cb-TDHFB which means they do not need 
Gram-Schmidt orthonormalization in the real time-evolution. We can expect the parallelization 
efficiency of the method is good, due to the relatively small communication-loss. 
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Figure 1. Neutron density distributions of XZ-plane in 20 O+ 20 O collision at i=(a)4.93, 
(6)98.66, (c)197.33, (d)295.99, (e)394.65 and (/)493.32 fm/c. The horizontal direction 
corresponds to Z direction. 
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Figure 2. Time-evolution of neutron gap-energy defined in Eq.([2D. Solid lines (red online) and 
dashed lines (blue online) indicate the gap energies of /-states in target and projectile nuclei, 
respectively. 



